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We consider a single molecule circuit embedded into solvent. The Born dielectric solvation model is 
combined with Keldysh nonequilibrium Green's functions to describe the electron transport proper- 
ties of the system. Depending on the dielectric constant, the solvent induces multiple nonequilibrium 
steady states with corresponding hysteresis in molecular current-voltage characteristics as well as 
negative differential resistance. We identify the physical range of solvent and molecular parame- 
ters where the effects are present. The position of the negative differential resistance peak can be 
controlled by the dielectric constant of the solvent. 



PACS numbers: 05.30.-d, 05.60.Gg, 72.10.Bg 

The use of molecules - either singly or in small ensem- 
bles - as the elements of electronic circuits holds sub- 
stantial promise in the fields of informational technology, 
biological and environmental nanosensors, and energy 
harvesting, [ij For the science of molecular electronics to 
be transformed into a technology it is not only important 
to fabricate stable molecular junctions but also to be able 
to efficiently control and manipulate their electric proper- 
ties. In the silicon-based microelectronic technology the 
gate voltage regulates the flow of electrons, but placing 
a third gate electrode has proven to be difficult in single 
molecular size devices. The negative differential resis- 
tance (NDR) also plays an important role in semiconduc- 
tor devices, because circuits with complicated functions 
can be implemented with significantly fewer components 
with its help. On the other hand, instead of copying the 
existing paradigms, such as, for example, gate voltage or 
resonant tunneling diode structure for NDR, the molecu- 
lar electronics create new and unique opportunities. The 
"wet" molecular electronics, where solvent controls the 
electric behavior of an electronic circuit, may open a new 
chapter in device engineering. Indeed, some molecular 
electronic devices already exploit the solvent around the 
molecule to modulate conductance through alteration of 
the charge state or polarizability of the molecule. 

Let us consider a "wet" molecular circuit - a molecule 
attached to two macroscopic metal electrodes and em- 
bedded into solvent (Fig. [T]). The total Hamiltonian is 



H = Hl + Hr + Hm + Ht + Hms- 



(1) 



The left and right electrodes contain free electrons and 
are described by the following Hamiltonians: 



(2) 



Here creates an electron with spin a in the single- 

particle state l/r oi the left/right electrode and ai^r/ra 
is the corresponding electron annihilation operator. The 



molecule is described by a single spin degenerate elec- 
tronic level with energy eq 
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The operator a\{aa) creates (destroys) an electron with 
spin a on the molecular level. The tunneling coupling 
between the molecule and electrodes is 



la 



ti{aj„aa + h.c) + 



+ h.c) 



(4) 



The interaction between the molecule and the surround- 
ing solvent, Hms, will be discussed below. We use nat- 
ural units in equations throughout the paper: h = ks = 
I — e| = 1, where — |e| is the electron charge. 

We describe the interaction between the molecule and 
the solvent based on the following simple model. The 
molecule is considered as a conducting sphere of radius 
R and the solvent is macroscopically uniform and char- 
acterized by dielectric constant e. The work needed to 
place charge qm on a conducting sphere in the dielec- 
tric environment is given bj^the Born expression for the 
dielectric solvation energy 
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where qs is the induced charge in the solvent {qm ~ 
—qs)- The model can be easily extended to the molecules 
of complex shapes (the so-called generalized Born model, 
which represents the molecule as a number of overlap- 
ping spheres of different radii). @ The (generalized) Born 
model is quite simple yet is very successful in comput- 
ing the electrostatic contribution to the solvation free 
energy. 01 The solvent dynamics is slow in comparison 
with the electron tunneling time scale. For example, the 
dielectric relaxation of the solvent is diffusive and occurs 
on the picosecond or slower time scales since dipolar sol- 
vent molecules generally respond to the change of the 
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FIG. 1. Schematic illustration of the model. The molecule is 
attached to two metal electrodes and surrounded by solvent. 
The solvent is described by uniform dielectric constant e. 



scales, so uj ^ 0.001 eV. For molecules interacting with 
the metal electrodes, F ~ 0.1 — 1 eV, which makes the 
static, effective mean-field (i.e., the static, average po- 
larization of the solvent) approximation (Eq. [6|) exactly 
valid for our case. 

Thus, in the Born approximation, the Hamiltonian 
Hm+Hi\is is exactly reduced to a spin degenerate single- 
level model with a local mean-field attractive interaction 
between electrons, which can be controlled by the dielec- 
tric constant of the environment. To describe electron 
transport through the system we use Keldysh nonequi- 
librium Green's-function formalism. [13, [Hj The exact 
nonequilibrium molecular population (N) and electric 
current J become: 



(N) = - f j^. ^L{u^)fLi^)+^Ri^)fRi^) 
^ ' TT J ^ {uj - e - 2A(w))2 
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(7) 



molecule junction charging state by rotating.]^ There- 
fore, we can assume that the induced charge qs corre- 
sponds to the average electronic population of the molec- 
ular junction. Then, the dielectric solvation energy can 
be directly associated with the interaction of the molecule 
with the surrounding solvent: 



Hms 



-U{e){N~6)i{N)~5), 



(6) 



where U{e) = 2^ (l ~ e) is an effective, local and sol- 
vent controlled electron-electron attraction, and N — 
^^aj-a^- The charge of the molecule due to nonequi- 
librium tunneling of electrons is {N — 5), while {S — (N)) 
is the corresponding induced charge in the solvent. The 
parameter S is the equilibrium molecular electronic pop- 
ulation which depends on the position of the molecu- 
lar level eo relative to the electrode Fermi energy £/. If 
Eq corresponds to the highest occupied molecular orbital 
(i.e., eo < then, without the applied voltage bias, the 
molecular level is double occupied and 6 = 2. If Eq is the 
lowest unoccupied molecular orbital (i.e. Eq > then 
the molecular level is empty in equilibrium and 6 = 0. 
It is known that such model Hamiltonians generally lead 
to bistable solutions. [1, 0] We emphasize that the model 
is not only applicable to the solvated molecular junction 
but also to the often employed experimental setting when 
the junction is embedded into isolating or semiconductor 
molecular film. In this case the surrounding molecular 
film can be considered as a macroscopic dielectric envi- 
ronment. 

The similar mean-field-type interaction between the 
molecule and the solvent (Eq. [S]) can be also obtained 
within the polaron model in the limit cj/F < < 1 (here w 
is the frequency of a characteristic vibrational mode cou- 
pled to the electrons and F is the broadening of the molec 



ular level due to coupling to the metal electrodes). [lO|,lll| 
In our case uj is related to the dielectric relaxation of the 
solvent, which occurs on the picosecond and slower time 
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Here /L/fl('^) = [1 + ef""^^-/^)/^]"^ is the Fermi-Dirac 
distribution for electrons in the electrodes, e = Eq — 
U{e){{N) — 5) is the effective energy of the molecular 
level, and A = -|- Aij, F = F^ -|- Fj,; are the real and 
imaginary parts of the electrode self-energy 
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AL/R{uj)-tTL/R{Lu). (9) 



The electrodes are modeled as a semi-infinite chain of 
atoms, characterized by the voltage-dependent on-site 



8 =-1.3eV 
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FIG. 2. Graphical solution of Eq. ((7|. The straight lines are 

given by equation (A'^) = — ^^jy -|- ^2 -|- jf^^ ■ Depending on e 

and eo there can exist one, three, or five fixed point nonequi- 
librium molecular populations. Filled circles represent sta- 
ble steady state populations, while open ones correspond to 
unstable fixed point solutions of Eq. ((Tj. Parameters: ap- 
plied voltage bias V = 1.5 eV, T = 300 K, Fq = 0.1 eV, 
R = 10 Bohr, e = 50 (dotted lines), and e = 3 (dashed lines). 
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energy tJLL,R = ±V/2 and the inter-site hoping parameter 
Vh = 2.5 eV. The expression for the electrode self-energy 



can be found, for example, in [ij]. The electrode band- 
is half filled, so the Fermi 



idth [jjL 



L/R- 



-2Vt 



h, ^J■L/R + ^Vh 



energy coincides with the one-site energy. The coupling 
between the left /right electrode edge and the molecule 
is taken to be VVftTo, where Fq = = ^r{^'r) is 

the maximal broadening of the molecular electronic level 
due to the coupling to the electrodes. Below we focus on 
the case when eo is lower than the electrode equilibrium 
Fermi energy. All our results also remain qualitatively 
valid when Eq is above the Fermi level. 

To compute the current, we first should determine the 
nonequilibrium molecular population {N). Since Eq. ([7]) 
is nonlinear, it generally has multiple solutions. Fig. [2] 
shows the graphical solution of this equation. As we see, 
likewise for the electron transport in the polaron model 
, depending on values of U (e) and the molecular level 
energy Eq, Eq. ([7]) can have one, three, or even five solu- 
tions (the nonequilibrium fixed points). These multiple 
solutions may or may not be steady states (i.e., the sta- 
ble fixed point). Following our method described in 15 1 
we obtain the stability matrix and analyze the real part 
of its spectrum to assess the asymptotic time behavior of 
the fixed points. We find that only the two outer and the 
middle solutions are stable in the five-solution case (right 
panel in Fig 12]); i.e., they correspond to physically real- 
izable nonequilibrium steady-state populations. In the 
case of three solutions (left panel in Figl2]), the middle 
solution is unstable and the other two fixed points are 
stable. We note that our approach is immune from the 
criticism that the observed multiple steady states are ar- 
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FIG. 4. Population-voltage and current-voltages characteris- 
tics. Parameters are T = 300 K, eo = — 1-6 eV, i? = 10 Bohr, 
e = 50, Fo = 0.1 eV. Three curves correspond to three possi- 
ble roots of Eq.((7]): solid line - upper root, dashed line - lower 
root, and dash-dotted line - middle root. 
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FIG. 3. The white domain corresponds to the values of 
£o and (7(e) where the multiple steady state solutions exist 
(Fo/y << 1). Dashed lines determine domain boundary and 
they are £o = -0.51/, eo = -U + 0.5^, £o = - 0.5^, 
eo = -2(7 + Q.W, and U = 1/27?. 



tifacts of the mean-field and electron self-interaction. [16j 
The effect of self-interaction is physically present in our 
case, since an electron in the molecule interacts with its 
own induced charge in the solvent. 

Let us now establish the range of key physical parame- 
ters - dielectric constant e, molecular size R, and molec- 
ular level energy eo , which allow the existence of multiple 
nonequilibrium steady states. For presentation purposes 
we assume that the molecular level broadening, Fq, as 
well as the temperature are much smaller than applied 
voltage V. Therefore the molecular population (Eq. [7]) 
(solid lines in Fig. [5]) can be approximated by a step like 
function of energy e. Then, we can readily determine an- 
alytically the conditions on Eq and U{e) when Eq. (O has 
only one solution. In Fig. |3] we show the domain where 
multiple steady states exist for the case Eq < 0. The case 
eo > can be considered in the same way, and the result- 
ing multistability domain is a mirror reflection of that in 
Fig. [3] across the abscissa axis. 

In Fig. Ill we show the behavior of the level population 
and the current as a function of applied voltage. Due to 
the presence of multiple steady states, both the popu- 
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FIG. 5. The NDR peak voltage value as a function of the 
dielectric constant for two different temperatures. Parameters 
are so = -2.0 eV, ii = 10 Bohr, To = 0.01 eV. Inset: NDR 
in current- voltage characteristic for e = oo. 



lation and the electron current demonstrate a hysteresis 
behavior. The width of the hysteresis loop is propor- 
tional to /7(e) and, therefore, it can be controlled by 
the dielectric constant. It should be emphasized that 
the solvent-induced hysteresis loop can be observed at 
moderate applied voltages where the molecular device is 
still mechanically stable. Moreover, the nonlinearity in 
the molecule-solvent interaction leads to NDR features 
in the current - voltage characteristic (the drop in the 
current represented by the dashed line at around 1 eV of 
applied voltage in Fig. S]). The NDR appears when one 
of the electrode chemical potentials crosses the position 
of the molecular level. Then, due to the subsequent shift 
in the level energy caused by the electronic population 
change, the level moves away from the current-carrying 
window between the chemical potentials. In the case of 
Eo < -0.5C/ (e) shown in Fig. H] the NDR takes place 
when we begin with the empty level. When Eq lays above 
— 0.5?7(e) (eg < 0) the NDR also takes place, but in this 
case we need to start from the initially fully occupied 
level. 

The NDR in the "wet" molecular circuit turns out to 
be sensitive to the dielectric constant of the environment. 
Figure [S] shows the dependence of the NDR peak position 
on the dielectric constant of the solvent. The increase of 
the solvent polarity shifts the peak toward the higher 
voltages. This effect is very robust. It does not require 
an artificial tuning of the model parameters and holds 
at very large ranges of temperatures. The temperature 
dependence of the NDR peak (inset in FiglS]) is consistent 
with experimental observations, and in contrast 

to the polaron model explanation of NDR Q does not 
require unphysical values for the parameters. 



We would like to comment here on the importance of 
the time scales. Depending on the relative time scales 
of measurements and transitions between stable fixed 
points, the multistability can result in merely noise as- 
sociated with the jumps between steady states or it can 
lead to hysteresis and NDR [l^ ■ To be experimentally re- 
solved the transition rate between multiple steady states 
should be smaller than the typical observation time. In 
our case the transition between steady states is deter- 
mined by the very slow diffusive reorganization of the 
solvent, which opens a possibility for experimental real- 
ization of the proposed effects. 

In conclusion, we have presented a theoretical model 
to describe the environmental control of the electron- 
transport properties of "wet" molecular junctions. The 
interaction between the molecule and solvent leads to ef- 
fective attraction between electrons which is governed by 
the dielectric constant of the surrounding solvent. The 
natural separation of electronic and solvent time scales 
makes the mean-field consideration exact for our model. 
We used Keldysh nonequilibrium Green's functions to ob- 
tain a nonlinear equation for molecular population and 
electric current. Depending on the dielectric constant, 
the inherent nonlinearity of molecule-solvent interactions 
induces multiple nonequilibrium steady states with cor- 
responding hysteresis in molecular I-V characteristics as 
well as NDR. We identify the physical range of solvent 
and molecular parameters which allows the appearance 
of multiple steady states. The temperature effects on the 
NDR peak are in qualitative agreement with the available 
experimental data. We demonstrated that the dielectric 
constant of the solvent can be used as a control parameter 
which regulates the position of the NDR peak. 
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